Bose-Einstein Condensation in Geometrically Deformed Tubes 
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We show that Bose-Einstein condensate can be created in 
quasi-one-dimensional systems in a purely geometrical way, 
namely by bending or other suitable deformation of a tube. 



It is well known [1] that there is no Bose-Einstein con- 
densate (BEC) of the continuous free Bose gas in one 
dimension as well as in tube-shaped domains DxR with 
a compact D C R 2 . On the other hand, there is an in- 
creasing experimental and theoretical interest [2] to the 
BEC in quasi-one-dimensional systems of trapped boson 
gases, in particular with the idea to understand prop- 
erties predicted by the Lieb-Liniger model [3] and the 
Girardeau- Tonks gas [4]. 

In this letter we want to draw attention to a way to 
create BEC in a quasi-one dimensional perfect Bose gas 
(PBG) based on geometrical properties of the recipient 
[5]. Dealing with a PBG we have naturally to describe 
first the one-particle spectrum. Consider first a straight 
round-tube-shaped recipient (cylinder) C(S r ) — S r x R 
of radius r > and the associated Dirichlet Laplacian [6] 
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The spectrum of t equals (Jj_i 00 )■> where 
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are the eigenvalues of the two-dimensional operator 
— (d% + dy) with the Dirichlet boundary conditions on 
the tube cross section S r . The lowest one is, of course, 
non-degenerate; for the circular cross section we have 

^ 1 = $m3o i r ~ 2 > where jo,i is the first zero of Jq(z) and 
the higher eigenvalues are similarly expressed through 
other Bessel functions zeros. The integral density of 
states corresponding to the spectrum of t is given by 



M(e) = 20rr(3/2) 
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where is Heaviside function. This gives for the PBG in 
the infinite tube the grand canonical total particle density 
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for temperature j3~ 1 > and chemical potential fj, < E\. 
Since the critical value p c {(3) '■= hm J u/'£ 1 p(/3,[x) of the 
density is infinite, there is no BEC of the PBG in this 
quasi-one-dimensional system. 



Several possibilities are known to make the critical 
density finite, provoking thus a BEC, by changing the 
one-particle spectrum. One is to replace the Dirichlet 
boundary condition by a "sticky" one, i.e. by a mixed 
condition with a positive outside gradient of the one par- 
ticle wave-function on the boundary [7]. One can also 
switch in an external local attractive potential producing 
bound state(s) below the inf o~(t) — E±, the continuum 
spectrum [8] . A less obvious way is a suppression of the 
density states at the bottom of the spectrum, leading to 
convergence of the integral (4) for p, = E\ , by embedding 
the PBG into a random external potential [9]. 

In the present note we are going to show that a geomet- 
rical deformation of the straight tube such as a simple 
local bending, even a gentle one, may produce the BEC 
of the PBG with the condensate localized in the vicinity 
of such a bend. We will also discuss other types of (lo- 
cal) geometrical deformations creating the BEC in these 
quasi-one dimensional systems as well the conditions un- 
der which the effect could be experimentally attainable. 

To this aim we have to recall the known results about 
one-particle spectra in deformed tubes. Consider first 
a local bending of the infinite cylinder C(S r ), so that 
the tube axis is a smooth curve which is straight out- 
side a compact region; we also suppose that the bent 
tube C*(S r ) does not intersect itself. If the deformation 
is nontrivial, C*(S r ) ^ C(S r ), it generates one or more 
eigenvalues [10,11] below the continuum threshold E\. 

Their number is finite and the lowest one of them is 
simple; the form of this discrete spectrum is determined 
the geometry of the deformed tube. Properties of such 
bound states in bent tubes are well understood: 
(al) If the tube is only slightly bent there is only bound 
state with energy and the gap E\ — t\ is proportional 
to ip 4 , where ip is the bending angle of C* , see [11]. 
(a2) The tube need not be straight outside a bounded 
region; it is sufficient that it is asymptotically straight in 
the sense that the tube axis curvature decays faster than 
|s| _1 where s is the arc length of the tube axis. 
(a3) The cross section need not be circular. In such a 
case, however, the shape is restricted by the so-called 
Tang condition imposed on the torsion [11]. It is satisfied, 
in particular, if C* = S* x [0, d], where S* is a bent planar 
strip [12]. 

(a4) The effect is robust. It also does not require the 
tube C* to be bent smoothly; the geometrically induced 
bound states exist also in sharply bent tubes [13]. 

Furthermore, bending is not the only way how to pro- 
duce geometrically a nonempty discrete spectrum: 
(bl) Another mechanism is a local change of the cross 
section. Protrusions and other deformations which do 
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not reduce the volume also give rise to an effective at- 
tractive interaction [14]. 

(b2) Similar effect comes from a tube branching. A cross- 
form tube is known to support a single bound state [15], 
numerous isolated eigenvalues arise in a skewed scissor- 
shaped cross with small enough angle [16]. 
(b3) One more example are two parallel tubes with a 
window in the common boundary, where the number of 
bound states is given by the window length [17]. 

The geometrically induced discrete spectrum is finite 
in asymptotically straight ducts. In the bent-tube case, 
for instance, it follows from the fact that the effective 
attractive potential falls off faster that |s|~ 2 , sec [11]. On 
the other hand, the actual form of the spectrum depends 
on tube geometry. In particular, if we deal with two 
or more well localized, identical, and mutually distant 
perturbations the eigenvalues cluster around those of a 
single deformation, with a split exponentially small w.r.t. 
the distance between the perturbations. This is proved 
mathematically for a tube with a pair of windows [18], 
however, various examples worked out numerically [19] 
suggest that such a behavior is generic. 

Consider thus such a geometrically induced discrete 
spectrum below the continuum threshold [20] in a de- 
formed cylinder which consists consists of one or several 
bound states with energies {e* < E\ : s = 1, 2, . . . } and 
cigenfunctions localized in the vicinity of the deforma- 
tion^). If one orders the bound-state energies naturally, 
ej < ei, < ... < Ei, then the domain of the allowed 
chemical potentials is given by the inequality p < e*. 
This implies that the critical particle density is bounded, 
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here Af*(e) is the integrated density of states for the de- 
formed cylinder. This opens a way to create the BEC. 

To make this effect transparent consider first a finite 
segment of the deformed cylinder C* L of length L. Then 
the Dirichlct Laplacian which plays here the role of the 
one-particle Hamiltonian, 



U{L) :-- 
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Dom(U(L)) = W^ 2 (C* L ), (6) 



has a purely discrete spectrum a(t*(L)) accumulating at 
infinity. Denote by Pi(t*(L)) the spectral projection of 
£*(£) for a Borel set I C R. It allows us to write the 
finite-volume integrated density of states corresponding 
to the operator (6) as 
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where \C* L \ is the volume of the segment C* L . By (7) the 
total particle density p* L ((3, p) in C* L acquires the form 
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Here {e*(L)} s >i and {e*j(L)}j>i are eigenvalues of the 
operator t*(L) divided into two groups with e* < e* for 
any s,j. The first one consists of those which converge to 
the eigenvalues of the infinite-tube operator t* as L —* oo; 
recall that all of them are monotonously decreasing with 
respect to L. The limit is naturally taken in such a 
way that the distance of the cut-offs from the deformed 
part(s) tend to infinity. On the other hand, {s*(L)}j>i 
are those eigenvalues which give in this limit the con- 
tinuous spectrum of t*. Combining this with the above 
stated properties of the one-particle spectrum we see that 
the first group {e*(L)} s >i is finite. Since \C* L \ — > oo as 
L — > oo, one gets then from (8) for the limiting density 
of states p*(f3,p) := liniL^oo Pl({3, p) the relation 
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provided /i < liniL^oo e\(L) = e\ < E\, uniformly in p, 
where N*{de) is the weak limit of the measure family 
{JVl(de)} L . In particular, the limit (9) implies 



p*(P,p)<p* c (P) :=/>*(/?, el) 



(10) 



for p < ej in view of (5). 

Now let the total particle density be p > p*((3). To 
show that this implies the BEC, let us analyze solutions 
{pl((3, p)} L of the equation following from (8), 



P = Pl(PiP) : 



(11) 



notice that they always exist and are bounded by e\{L) 
from above. Using (8) and (9) one gets for p < p*((3) 



lim p L (f3,p) = p(f3,p) < . 

L — >oo 



(12) 



On the other hand, for p > p*{(3) we rewrite eq. (11) as 

-— V 1 

9 \r*\ 



\C*\ ^ J3( S '(L)-n L (l3,p)) 

1 Ll {*;w} e 



1 



- y 

r* I ^ 



1 



\ e f3(e*(L)-p L (0,p)) _ l 



(13) 



(i) Suppose that the bound-state energies {Cs(-k)} s >i ver " 
ify the following conditions as L — > oo, 



Ei - 8 > e[{L) for some 5 > 



and 



e:(L)-eI(L)>o|Ci 



,a-l 



a > 0, a > 0. 



This is true, in particular, for a tube with a single bend 
or protrusion and two cut-offs moving away of it, where 
in view of the norm-resolvent convergence of — Aj-f to 
— A^,* the eigenvalue e\(L) tends to e\ < E\ and the 
eigenvalue difference to a nonzero limit. 



Since Pl{P,p) < e*{L), we can use the uniform con- 
vergence (9) to obtain the asymptotics of the solution of 
cq. (13) as L — > oo, namely 
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0(p-p*M) IQ 



o (\c* L 



(14) 



This means that hin^^oo Pl((3, p) = e*, and since \C* L \ = 
O(L) as L — > oo, the thermodynamic limit in (13) gives 
rise to BEC at the lowest level only 
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The same is true for other tube geometries, for instance, 
for a tube with several well distinguished bends or "bub- 
bles" , as long as the limit L — > oo means that cut-offs 
move away in the asymptotically straight parts. 

(ii) The situation may change when the thermody- 
namic limit L — > oo is more complicated and involves a 
local change of the geometry as well. As a model exam- 
ple, consider again a tube with a finite number n > 1 of 
well distinguished identical bends, but suppose now that 
the distances between them increase also with increas- 
ing L. As we have recalled above, the first n eigenvalues 
{e*(L)}™ =1 then cluster being exponentially close to each 
other with respect to the separation parameter L, 



){L)-e\{L)<Ce 
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1 < s <n. 



(16) 



for some positive C, a. By (16) the asymptotics of the 
solution to eq. (13) is again of the form (14). This means 
that the limit in (13) gives rise to BEC equally frag- 
mented into the group of n levels {e^{L)}™ =1 which are 
almost degenerate being exponentially separated, 
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for 1 < s < n. This fragmentation is called a type-I 
generalized BEC , contrasting with the case of the infinite 
fragmentation known as the type- II generalized BEC, see 
[1], [21] and also [22]. What is important is that this 
condensate is separated from the continuum spectrum 
by a finite energy gap which makes it more stable than 
the conventional BEC, see the discussion in [8]. 

To analyze localization properties of the geometrically 
induced BEC we employ the PBG one-body reduced den- 
sity matrix with the kernel Pl(P, P', x > V) given by 

r°° i 

\c* L \ u* L {de) r e , L {xW StL {y) > ( 18 ) 

J —oo ° * 

where {ipt l} are the normalized eigenf unctions of the 
operator (6). The diagonal part of the matrix (18) is the 
local particle density, since by (8) 



c; 



dxp L (f3,p;x,x) = \Cl\p L (P,n) 



(19) 



is the total number of particles in C* L . In fact a rele- 
vant quantity for the BEC space localization is the local 
particle density per unit volume pl{P, p; x) given by 



r°° i 

J —oo c 



(20) 



Indeed, since in the limit only the cigenfunctions cor- 
responding the eigenfunction family {"0s} s >i C L 2 (C*) 
related to the infinite-tube bound states is preserved, in 
contrast to all the others which are extended, we get 

f°° 1 

p{(3,p;x):=Y^ J M* L {de) _ ^ ^(^(i) 

= (p-p*(/3))i^)€ 1 (x) (21) 

for p > p*((3) in the case (15). For the fragmented BEC 
(17) one gets for p > p*(/3) similarly 



p(f3,p;x) 



(22) 



S = l 



it is obvious that p(p, p;x) =0 holds when p < p*([3). 

Thus in contrast to the BEC in the translation- 
invariant case, the L 2 -localized condensation in the 
bounded states correspond to an infinite accumulation of 
the local particle density defined by (18). On the other 
hand, by (15) and (17) in combination with (21), (22), 
the total density of particles condensed in the bounded 
states, 



P~ P*c(P) = J c dxp((3,p; : 



(23) 



is finite. 

After this theoretical analysis let us ask about chances 
to observe the described type of the BE condensation in 
an experiment. We are not going to discuss technically 
the ways in which the Bose gas can be confined in a 
geometrically deformed tube, generally we have in mind 
either a modification of the existing elongated traps [23] 
or using hollow optical fibers [24] as suggested in [25]. 

The important parameter is the tube radius which for 
both the elongated traps and hollow fibers can be made as 
small as r « 5 pm which determine the threshold energy 
by the mentioned expression E\ = jo.i r ~ 2 1 where M 
is atom of the mass in question. Let us further intro- 
duce the relative gap size by 7 ro i := (Ei — e\)/Ei. This 
quantity can theoretically reach the value about 0.39 in a 
smoothly bent tube [26] , but a typical value for a bending 
angle of 90° and more is around 7 rc i « 10 _1 , see [19]. 

Two conditions must be satisfied. First of all, the 
bending-induced gap 7 re i-Ei must be much larger than 
the effect of the finite length L of the recipient. The lat- 
ter is characterized by the first longitudinal eigenvalue 
(tt/L) 2 . It is clear that the trap must sufficiently elon- 
gated to fulfill the condition, roughly speaking L/r > 20 
is sufficient. This is true in both situations mentioned 
above, for hollow fibers it can be even much better. 

More complicated is the question of thermal stabil- 
ity; the said gap must be larger than the energy k^T. 
This determines a critical temperature above which the 
bending-induced BEC is likely to be destroyed by ther- 
mal fluctuations. Using the value r«5 pm we find that 



T crit » 5.4 x 1(T 8 ^ , (24) 

where Z is the atomic number in question. Hence lighter 
nuclei are preferable; there is one order of magnitude dif- 
ference between Li 7 and Ru 87 . With above estimate in 
mind, however, even for the light ones it is difficult to per- 
form the measurement in the available nanokelvin condi- 
tions. On the other hand, the effect is not too far from 
the experimental reach; recall for instance that squeezing 
the transverse size to a one micron radius would enhance 
the critical temperature (24) by the factor of 25. 

Let us finally mention one more feature of such a ge- 
ometrically induced BE condensate. The ground-state 
wave function into which will the atoms of the Bose gas 
condensate - cf. (21) - is exponentially localized away 
of the bend; for examples of such wave functions see [f 9] 
and the literature mentioned there. This means that, 
e.g., bending a tube would mean not only the conden- 
sation but also that the condensate will squeezed in the 
bend the more the larger the bending angle is. 

In conclusion, we have demonstrated here a purely geo- 
metric way to achieve a stable BEC based on local defor- 
mations of a tube-shaped recipient and we discussed ways 
in which the effect could be experimentally observed. 
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